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COMPUTERIZED ARRHENIUS 
RELIABILITY EXTRAPOLATION TECHNIQUES 


By Charles R. Toye 
Electronics Research Center 


SUMMARY 

This paper presents computerized mathematical techniques and algorithms 
that are useful in obtaining extrapolated component failure rates derived from the 
Arrhenius equation. The mathematics of the Arrhenius equation is first reviewed. 
A computer program for the calculation of an important constant in the equation 
is then given, and finally, a computerized extrapolation procedure is illustrated. 

INTRODUCTION 

The Arrhenius relation plays a significant role in the theory of physical- 
chemical reaction kinetics and has been used extensively for the prediction of 
chemical reaction rates as a function of temperature. In the semiconductor 
industry it is the most commonly used statistical technique for the extrapolation 
of failure rates. To the extent that semiconductor failures are traceable to 
temperature-generated, physical -chemical reactions of the device materials 
involved, the Arrhenius relation provides an adequate basis for failure rate 
extrapolation . 

In its exponential form the relation is given by 

-C/T 

F 2 = F l e <>> 

whe re 

Fj = the failure rate (expressed as a decimal per 1000 hours) 
at maximum test temperature T^ 

= the failure rate (similarly expressed) at some lower 
temperature T^ 

C - the Arrhenius accele ration factor related to the particular 
device considered, measured in degrees Kelvin 

T ~T , ‘ " 

1_1 1 12 _ 

= = — = — =- .where T and T 'are .measure in-degrees Kelvin. 
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By taking common logarithms of Eq. (1), it can be shown that the logarithm 
of the unknown failure rate can be expressed as a linear function of 1/T. Thus, 
we have : 




= log 


Now: 


and 


log F 2 = log ( F 1 e' C/T ) 


log F ^ = log Fj + log e 


■C/T . 


, -C/T _ , -C/T 
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log F 2 = log Fj - Q ■ 


2. 3-03 


( 2 ) 


A sample problem is now given to illustrate its application. 

PROBLEM : Suppose 100 units of a particular device are tested for 1000 hours at 

300 C and a failure rate of 10%/ 1000 hours is observed. A second sample is 
tested at 125 C with an observed failure rate of 0.01%/ 1000 hours. Determine 
the failure rate (%/ 1000 hours) at 60°C. 


SOLUTION: To use the Arrhenius relation a value for C (Arrhenius acceleration 

factor) must first be determined. For this, the logarithmic form of the relation 
given by Eq. (2) is best suited: 


log F 2 
-C 
C 
C 


Substituting known values: 


F =0.1 

Tj = 300°C = 573°K 

f 2 = 0.0001 

T 2 = 1 25°C = 398°K 


* In changing from Centigrade to Kelvin scales it is necessary to add 27 3 
to the Centigrade reading. 


= log F. 


C 


1 2.303 

= 2. 303 T 


(t) 


(log F 2 - log Fjj 
= 2. 303 T ^log Fj - log F^ 


= 2. 303 


(^S)( log F ‘ ' log Fz ) ' 
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we have: 


C = 2. 303 
C = 2.303 


(HH-iH) ( log 1 - 108 - 0001 ) 

/573 xJ98\ / , +4 ) 

\573 - 398/ \ / 


C ~ 9003. 57°K 

Once a value for C is obtained, the extrapolated failure fate at 60°C 
can be determined by either of the following two methods. Let F = 0. 1 and 
T = 573°K, with C = 9003. 57°K and T z = 60°C = 330°K. 


METHOD 1 

Using Eq. (1) and substituting known values we have: 



F 2 = 0.1 (e- 10 )(e- 132471 ) 


F 2 « 0. 1 (.0000454) (.265883) 

F 2 w 0.0000012071 . 00012%/ 1000 hours. 


METHOD 2 


U'sing Eq. (2) and substituting known values we have: 


log F 2 

1 IT C 

- log F. - 

1 2. 303 

log F 2 

, , 9003.57 

= l0g ‘ 1 - 2.303 

log F 2 

j 900 3.57 / 

2. 303 \ 

log F 2 

K -1 -4.9174 = -5. 

F 2 

Rj 1. 209 x 10“ 6 = 0. 


/ T l- T 2 \ 

V T 1 T 2 / 


/ 573 - 333 
\^573 x 333 

573 - 333\ 
573 x 333 ) 


0. 00012%/ 1000 hours. 
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FACTOR PROGRAM 


The Arrhenius relation in its logarithmic form (Eq. (2)) is graphically 
presentable on semi-logarithmic paper as a linear curve. Linear curve fitting 
techniques, e.g. , the "Least Squares Method, " would require data derived from 
several test temperature conditions. However, because of the high cost in the 
semiconductor industry, the experimental design is generally limited to two 
test conditions. 

Contained herein (Figure 1) is a computer program (FACTOR) using the 
IBM 7090/7094 IBSYS Operating System, FORTRAN IV language*, it is designed 
to calculate the value of the Arrhenius acceleration factor C necessary for the 
extrapolation of failure rates. In essence, the program utilizes the same 
method used in the sample problem above to calculate C. Given any two test 
temperatures ( C) and their respective failure rates (%/ 1000 hours), or the 
information necessary to determine these failure rates, i.e. , number of 
failures, sample size, and length of test (hours), the program used the Arrhenius 
relation (Eq. (2)) and the computer's logarithm subroutine to calculate the 
Arrhenius acceleration factor C ( K) for the particular device under test. 

The program is set up to accept test data in either of two forms. Form I 
consists of two test temperature ( C) and their respective failure rates 
(%/1000 hours). Form II consists of two test temperatures ( C) and their respec- 
tive number of failures, sample sizes, and hours tested. 

INPUT SETUP AND PROCEDURE 

All input data must be keypunched on tabulating cards. The data deck 
consists of at least two data cards for each device for which an acceleration 
factor is to be calculated. For each device, the first card(s) contains actual test 
data; the last card, descriptive information, i. e. , the kind of device and its type 
numbe r . 


To input test data of the kind in either Form I or Form II, the first data 
card(s) for each device considered, except the last, must be either of the form 

Column 23 

$SPEC TEMPI =300 , TEMP2 =200, FR1 =10, FR2 = 7$ 


or 

Column 23 

$SPEC TEMP 1= 300 , TEMP2=200, REJ 1= 10 , REJ2-7, SIZE1= 100 , SIZE2=1000 , 
HRS1 =1000, HRS2 =1000$ 

depending .-n whether test data is of Form I or Form II, respectively. The mean- 
ing of the variables is as follows: TEMPI is the higher of the two test tempera- 

tures (°C); TEMP2 is the lower of the two; FR1, REJ1, SIZE1 and HRS1 are, 
respectively, the observed failure rate (%/ 1000 hours), number of failures, 
sample size and hours tested at TEMPI; FR2, REJ2, SIZE2, and HRS2 are 
similarly defined for TEMP2. 

In setting up the "SPEC" data card(s), column I is left blank; a $ must 
appear in column 2, immediately followed by the word "SPEC", immediately 
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followed by one or more blank characters. The variable names are then assigned 
constant values, with all data items separated by commas. Use of a comma 
following the last item is optional. If more than one card is needed for "SPEC" 
data, data items begin in column 2 and the last item of each card, except the last, 
must be a constant followed by a comma. The end of a group of "SPEC" data is 
signaled by a $ either in the same card as the word "SPEC", or anywhere in any 
succeeding card except in the column 1 position. 

The descriptive data card for each device considered is of the form 
Column 1 18 

TRANSISTOR 2N2651 

Columns 1-12 are rese rved for the kind of device, e.g., transis tor , diode , etc . , 
entered left -jus tified ; columns 13-18 are reserved for the type number, e.g. , 
2N2651, 1N696, etc. , entered right-justified. 

Calculation of an acceleration factor for more than one device is accom- 
plished by simply supplying additional "SPEC" and related descriptive data cards. 
Input data for the last device for which an acceleration factor is to be calculated 
must be prepared in either of the following forms: 

C olumn 12 3 18 

$SPE C TEMP 1 = 300, TEMP2=200, FR1=10, FR2=. 7, FINAL=1$ 
TRANSISTOR 2N2651 


or 


Column 123 18 

$SPEC TEMP 1 = 300, TEMP2=200 , R E J 1- 1 0 , REJ 2=7, SIZE1=100, 

SIZE2 =1000 , HRS 1 =1000, HRS2=1000, FINAL=1$ 

TRANSISTOR 2N2651 

again depending on whether test data is of Form I or Form II, respectively. Note 
that the input setup for the last device differs from that for previous devices only 
in the addition of "FINAL=1" to the "SPEC" card(s); this causes execution to 
te rminate. 


Output is the form: 

DEVICE TYPE ACCELERATION FACTOR 


TRANSISTOR 2N345 
TRANSISTOR 2N2651 


9234. 71369288 
7208. 66473389 


The following is a listing of the program proper (Figure 1) followed by 
sample input (Figure 2) and output (Figure 3) listings and a sample input card 
setup (Figure 4). 
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FIGURE 1 


ACCELERATION FACTOR (FACTOR) 
PROGRAM LISTING 


SJOB CABRAL-NASA 

SEXECUTE IBJOB 

S IB JOB 

SIBFTC FACTOR 

NAMELI 5T /SPEC /TEMPI ,TEMP2 »FR1 *FR2»RE Jl »RE J2 » S I ZE 1 »S I ZE2 * HR SI ,HRS2» 
IFINAL 
WRITE( 6,5) 

5 FORMAT ( 1H1 , 30X , 6HDE V I CE , 1 OX , 4HT YPE , 2 OX , 1 9HACCELER AT I ON FACTOR) 

WR ITE( 6, 12 ) 

12 FORMAT ( 1HO ) 

9 TEMP1=0 
TEMP2= 0 
FR 1 = 0 
FR2 = 0 
RE J1=0 
RE J2=0 
S I ZE 1 = 0 
S I ZE2= 0 
HRS1=0 
HRS2=0 
F I NAL = 0 
RE AD ( 5 ,SPEC) 

READ ( 5 , 2 ) DEVI ,DEV2,TYPE 

2 FORMAT (3A6) 

I F ( FR 1 )3,A,3 

3 FR1=FR 1/100, 

FR2=FR 2/100. 

A FR1=REJ1/(SIZE1*HRS1 1*10. **5 
FR2 = REJ2/< SlZE2*HRs2)*10.**5 
X=TEMPl+273. 

Y=TEMP2+273. 

T = ( X— Y ) / ( X*Y ) 

W=( (ALOGIO (FR2 ) -ALOGIO (FR1 ) >*2.303) /T 
C=-l .*W 

WR I T E ( 6,6) DEVI »DEV2 ,TYPE,C 

6 FORMAT (30X,A6,A6,AX,A6,18X,F18.8) 

I F ( F I N AL-1 • ) 9 , 8 , 9 

8 STOP 
END 


FIGURE 2. 


SAMPLE INPUT LISTING 


$ DAT A 

5 SPEC TEMPI =300,TEMP2=2O0,FRl=10,FR2=.7$ 

TRANSISTOR 2N2651 

SSPEC TEMP1=300,TEMP2=200»REJ1=10»REJ2=7 , S I ZE 1 = 100 ,S I ZE2=lOOO »HRS 1 = 1 000 , 
HRS2=1000,FINAL=lS 
TRANSISTOR 2N2651 
$ STOP 


FIGURE 3. 


SAMPLE OUTPUT LISTING 


DEVICE TYPE ACCELERATION FACTOR 

TRANSISTOR 2N2651 7208 .66A73389 

TRANSISTOR 2N2651 7208.66A73389 
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The two acceleration factors outputted (Figure 3) are identical and are, 
in fact, for the same device. Checking the input data (Figure 2), one finds 
that the first of the two identical lines of output corresponds to input data of 
the kind described by Form I, while the other corresponds to that described by 
Form II. Data inputted the second way required the computer to first calculate 
failure rates (which, in this case, matched exactly those failure rates just 
previously inputted directly) before determining an acceleration factor. The 
duplication of output lines was purposely arranged in the data deck to illustrate 
the program's duality of input modes. 

EXTRAP PROGRAM 

The purpose of the above FACTOR program (Figure 1) is to calculate the 
value of the Arrhenius acceleration factor C, which in turn can be used in the 
extrapolation of failure rates for the device under test. A program (EXTRAP) 
has been written to perform this latter task. As with FACTOR, EXTRAP uses 
the IBM 7090/7094 IBSYS Operating System, FORTRAN IV language. 

Specifically, given data based on performance under high temperature 
testing, i. e. , temperature ( C) , sample size, length of test hours, number of 
failures and acceleration factor ( K) , the EXTRAP program will output an 
extrapolation table(s) which lists temperatures ( C) in descending order, begin- 
ning with the given high test temperature and decreasing by decrements of 5 C. 
The program is written to accommodate any high test temperature not greater 
than 200 C (the temperature at which low- and medium-power silicon trans- 
istors are perhaps most frequently life-tested). For each temperature between, 
and including this maximum test temperature and a minimum temperature not 
lower than 25°C, the table includes a failure rate (%/ 1000 hours) for each of 
five confidence limits of 10, 60, 70, 80 and 90 per cent. (These limits were 
calculated by means of the "Chi-square distribution".) Thus, on the basis of 
high temperature performance test data, the reliability of components under 
"in-use" conditions, through the Arrhenius extrapolation procedure contained 
within this program, can be predicted and presented in tabular form. 

INPUT SETUP AND PROCEDURE 

Input data is keypunched on tabulating cards and is set up as outlined 

below. 


The first nine cards are reserved for general descriptive and/or identify- 
ing information, punched in columns 1-72. Data presented on these first nine 
cards will appear alone on the first page of output and serve as a title page for 
the extrapolation table(s) to follow. 

The next card contains the actual test data that forms the basis for the 
extrapolation of failure rates. The same general procedure is followed that has 
been described for setting up the "SPEC" cards in the preceding FACTOR pro- 
gram, i. e. , starting with a $ in column 2, etc. Test data is thus presented in 
the following form: 

Column 23 

$SPEC SAMPLE=332, HRS=1000, TEMP=200 , NF =1 1 , C=7500, N=6$ . 
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The meaning of the variables is as follows: SAMPLE refers to the sample size 

or number of units of the device tested, HRS refers to the length of the test 
{hours), TEMP refers to the temperature (°C) in the range 25 C to 200°C at 
which the test was conducted, NF is the number of failures observed, C is the 
acceleration factor (°K) and N is simply a number, the value of which controls 
the continuation or termination of execution. When N is less than or greater 
than zero, execution continues; when N equals zero, execution is terminated. 

The next three cards are reserved for descriptive information, punched 
in columns 1-72. These data will supply a heading for the extrapolation table 
based on the immediately preceding "SPEC" card. 

The EXTRAP program can accommodate more than one set of test data 
and related descriptive information and can thus output more than one extra- 
polation table on a single run. If this is desired, one need only supply the 
additional "SPEC" cards each followed by three descriptive cards. The last 
"SPEC" card must assign to N a value of zero (see Figures 5 and 6). 

The following is a sample input listing (Figure 5), its corre sponding 
card setup (Figure 6), and its output (Figure 7). The tables in Figures 7b and 7c 
correspond to sampling plans in Table C-I of the MIL-S-19500, indexed accord- 
ing to the lambda value, sample size, and number of failures permitted (accept- 
ance number). In Figure 7b, lambda equals 10, the sample size is 178, and the 
acceptance number is 12. In Figure 7c, lambda equals 5, the sample size is 
332, and the acceptance number is 11. 

In Figures 5 and 6, cards 13 and 17 are blank because only two lines of 
descriptive information were sufficient to identify each table and three descrip- 
tive cards are required by the program as it is written. In output (Figures 7b 
and 7c), the blank card causes a space between the two lines which head each 
table. 


If it is so desired, one may use blank cards similarly in the first set 
of nine descriptive cards used as a title page for the tables, as long as the 
number of cards in this set (blank and/or punched) totals nine. 

The Arrhenius equation was developed by S. Arrhenius in 1889 as an expression 
for the influence of temperature on reaction velocity. A derivation of this relation- 
ship may be found in Reference 1 . 

Variations of the Arrhenius relation for use in semiconductor reliability may 
be found in Reference 2 . 

A brief discussion of the Arrhenius equation and its limitations may be found in 
Reference 3 
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FIGURE 5. 

SAMPLE INPUT LISTING 


SDAT A 

RELIABILITY EXTRAPOLAT ION BASED ON ARRHENIUS RELATION 
SEMICONDUCTOR DEVICES 

FAILURE RATE VARIATION WITH TEMPERATURE 

EXTRAPOLATION - 5 CONFIDENCE LIMITS (10* 60, 70, 80» 90 PERCENT) 
INFORMATION KNOWN - NUMBER OF FAILURES, MAXIMUM TEMPERATURE, 

SAMPLE SIZE, LENGTH OF TEST (HRS.), ACCELERATION FACTOR 
INFORMATION LEARNED - 5 FAILURE RATE CONFIDENCE LIMITS F.OR EACH TEMPERA- 
TURE WITH TEMPERATURES DECREASING BY DECREMENTS OF 5 DEGRtES 
NASA-ERC (CQS) DATE - AUGUST 1, 1967 

$SPEC SAMPLE=178,HRS=1000 .TEMP* 200 »NF=12»C=7500,N=6$ 

LAMBDA = 10 

f* - blank card T 

MIL-S-19500 SAMPLING PLAN - TABLE C-I 
$ SPEC SAMPLE=332»HRS=1000»TEMP=200»NF=ll»C = 7500,N = 0$ 

LAMBDA = 5 

r blank card | 

MIL-S-19500 SAMPLING PLAN - TABLE C-l 
S STOP 
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FIGURE 6 
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c» HRS=10Q0» TEMP=ROO, HF=1 1 j C=7500» N=0$ 


WEammmamuaMaasaau 


SPEC SAMPLE=17A*HRS=in00,TEMP=P00,NF=lR»C=7500»N 


ASA-ERC (CQS) DA jr - AUGUST t, 1% 


/ TURF i-UTH TEMPERATURES DECREES I HG BY DECREMENTS OF 5 DEGREES 


NFORMATION LEARNED - 5 FAILURE RATE CONFIDENCE LIMITS FOR EACH 


/ SAMPLE SIZE* LENGTH OP TEST <HRS.>, ACCELERATION FACTO 


INFORMATION KNOUN - NUMBER OF FA I LURES j MAXIMUM TEMPERATURE* 


.XTRAPOLAT i ON - 5 CONFIDENCE LIMITS <10* AO* 70, 80, 00 PERCENT) 


A I LURE RATE VARIATION UITH TEMPERATURE N 
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FIGURE 7. 
SAMPLE OUTPUT 


RELIABILITY EXTRAPOLATION BASED ON ARRHENIUS RELATION 
SEMICONDUCTOR DFVTCES 

FATLURF PATE VARIATION WITH TEMPERATURE 

EXTRAPOLATION - 5 CONFIDENCE LIMITS (10* -60* 70* 80* 90 PERCENT) 
INFORMATION KNOWN - NUMRFR OF FAILURES* MAXIMUM TEMPERATURE* 

SAMPLE SIZE* LENGTH OF TEST (HRS,)* ACCELERATION FACTOR 
INFORMATION LEARNED - 5 FAILURE RATE CONFIDENCE LIMITS FOR EACH TEMPERA- 
TURE WITH TEMPERATURES DECREASING 8Y DECREMENTS OF 5 DEGREES 
4MASA-ERC (COS) DATE - AUGUST 1* 1967 


-Figure 7a. 


LAMBDA = 10 


M IL-S-195Q0 SAMPLING PLAN - TABLE C-I 


f r NG T H °P rrrcj 
1 Of'"', 


SA M P| F SI 7 F 
1 7R. 


Nfl. PF FAILURES 
12 


MAXIMUM TEMPERATURE 

200 . 


ACCELERATION FACTOR 
7500. 


RELIABILITY EXTRAPOLATION - FAILURE RATE (PERCENT/1000 HOURS) 


IIP F 

rFMTTPo ^r)P 

confiofncf 

lo 

LIMITS (PERCENT) 
60 

70 

80 

90 

*>00. 

4.859550 

7.640449 

8.202247 

8.932584 

1C.OOCOOO 

lo<5. 

4.102275 

6.449820 
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